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ABSTRACT 

We present observations spanning 355 orbital phases of HAT-P-7 observed by Kepler from May 
2009 to March 2011 (Ql-9). We find a shallower secondary eclipse depth than initially announced, 
consistent with a low optical albedo and detection of nearly exclusively thermal emission, without 
a reflected light component. We find an approximately 10 ppm perturbation to the average transit 
light curve near phase -0.02 that we attribute to a temperature decrease on the surface of the star, 
phased to the orbit of the planet. This cooler spot is consistent with planet-induced gravity darkening, 
slightly lagging the sub-planet position due to the finite response time of the stellar atmosphere. The 
brightness temperature of HAT-P-7b in the Kepler bandpass is Tb = 2733 ± 21 K and the amplitude 
of the deviation in stellar surface temperature due to gravity darkening is approximately —0.18 K. 
The detection of the spot is not statistically unequivocal due its small amplitude, though additional 
Kepler observations should be able to verify the astrophysical nature of the anomaly. 
Subject headings: Planets and satellites: individual: HAT-P-7b — Planets and satellites: fundamental 
parameters — Stars: Individual: HAT-P-7 — Stars: atmospheres — Eclipses 



1. INTRODUCTION 

Since launching in 2009, the Kepler Space Telescope 
has provided unparalleled photometry of more than 
160,000 stars, with the primary goal of detecting tran- 
siting extrasolar planets. By providing continuous short- 
cadence photometry across a single field of view, Ke- 
pler has succeeded in detecting thousands of new c an- 
didate transiting p lanets (see IBorucki et al.l 120111 or 
iBatalha et alj 120121 for a recent review) , but it has also 
provided extremely high-precision photometric time se- 
ries for three extrasolar planets that were known prior 
to launch: TrES-2b, HAT-P-llb and the hottest of the 
bright Kepler planets, HAT-P-7b. These targets, due 
to the brightness of the parent stars and their well- 
constrained orbital parameters, represent an opportunity 
to search for transit-curve features that are unobservable 
for any other objects at present, and which can inform 
us about the characteristics of the planet-star system. In 
this Letter we report on our analysis of HAT-P-7b, the 
brighter and hotter of the three objects. 

HAT-P-7b was originally discovered by iPal et all 
(2008) in the HATNet project. This transiting hot 
Jupiter with mass 1.7 Mj orbits an F6 star with an or- 
bital semi-major axis only four times greater than the 
radius of its host star. HAT-P-7b is one of the hottest 
known transiting exop lanets in the Kepler fie l d, an d 
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atmospheric models by iMadhusudhan fc Seagerl 
iSpiegel fc Burrows! (|2010f ) and lChristiansen et ail 
hereafter denoted C10, suggest strong evidence for ther- 
mal inversion. 

As a known bright exoplanet target, HAT-P-7 was 
given high priority by Kepler and has been observed 
using a one-minute cadence nearly continuously since 
May 2009. Using only the first 10 days of observations, 
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IBorucki et all (|2009l ) were able to detect the secondary 
eclipse at optical wavelengths, finding an eclipse depth 
of 130 ±11 parts per million ( ppm) as well as orbital 
phase variations. iWelsh et al.l (|2010f) determined that 
these phase variations were dominated by perturbations 
due to ellipsoidal variations in the shape of the host star 
caused by the gravity of the planet, resulting in brighten- 
ing of the system on either side of the secondary eclipse 
with an amplitude on the order of ^60 ppm (Jackson et 
al. 2012, hereafter denoted J12). 

Each of these studies was limited by the amount of 
data available at the time - the most recent study by 
J12 analyzed data up through Quarter 2. With the re- 
cent release of newly re-calibrated data spanning Q1-Q9, 
an updated analysis of the full set of transits and eclipses 
(326 transits and 355 eclipses in all, after selection cri- 
teria are applied) is warranted. In this work we focus 
on analyzing both the variability of the eclipse depths, 
as well as the characteristics of the mean transit and 
eclipse curves. These measurements allow us to place 
improved constraints on the presence of cloud variability 
and planetary albedo, and to search for additional effects 
of the planet on the host star. In |J5]we describe our ob- 
servational data set, in Sj3]we describe our transit fitting 
methods and analysis methodology, and in S}4] and $5] we 
discuss and summarize our results. 

2. OBSERVATIONS 

Short-cadence (^60 s) photometry of HAT-P-7 was 
recorded by Kepler from 2009 May 13 - 2011 June 26 or 
Quarters 1-9. More than one million photometric mea- 
surements produced by the Kepler data analysis pipeline 
were retrieved. We examined fluxes calculated from both 
SAP (Simple Aperture Photometry) and PDCSAP (Pre- 
search Data Conditioned Simple Aperture Photometry). 
The conditioning process, a standard component of the 
Kepler data pipeline, seeks to remove systematic effects 
introduced by the sp acecraft while retain ing interesting 
astronomical signals (| Jenkins et al.l [20~ioh . The results 
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Fig. 1. — Upper: Composite transit light curve of HAT-P-7b 
from 326 transits. The >400,000 photometric measurements in 
the range of phase from -0.1 to 0.1 are averaged in 1-minute (334 
photometric measurement) bins. Error-bars are omitted for clar- 
ity. M iddle: Residuals after removing the best-fit Mandcl & Aeol 
(2002) light curve model using Levenburg-Marquardt least squares. 
The c urve is the best - fit m odel light curve for the anomaly 
in the Mandcl & Agol (2002) model caused by the planetary 
transit across the gravity-darkened spot on the surface of the 
star. Lower: Residuals after removing the minimum-x 2 gravity- 
dar kened spot model (circ les) and the expectation after removal of 
the Mandcl & Agol (2002) and gravity-darkened spot models (line 
at ppm). The standard deviation of these final residuals is 11.6 
ppm. 

we show here were calculated with the PDCSAP fluxes, 
since we expect these data to contain less uncorrected 
systematic effects, and thus to yield more accurate re- 
sults. We carried out the analysis with both types of 
flux and found that the results were consistent with each 
other, despite the additional corrections imposed on the 
PDCSAP data set. 

Each successive orbital phase of observations of the 
planet was considered for inclusion in our analysis in- 
dividually. We rejected orbital phases that contained 
uncorrected systematic effects after the Kepler pipeline 
removal and any observation that did not record a com- 
plete time series on phase -0.2 to 0.2 (near the transit) or 
0.3 to 0.7 (near the eclipse), where the phase is defined 
on 0-1, with zero centered on mid-transit. 

3. ANALYSIS 
3.1. Primary Transit 

Using the orbital parameters determined by C10, we 
calculated an initial estimate for the mid-transit time 
for each transit. Model light curves of the form de- 
rived by iMandel fc Agol |2002) were then fit to each 
transit over orbital phase -0.2 to 0.2 using Levenburg- 
Marquardt least-squares. We fit for the phase angle of 
mid-transit, ratio of planetary to stellar radii (Rp/R*), 
ratio of the semi-major axis of the planet to the radius 
of the star (a/i?»), inclination (i), zero-point flux and 
linear and quadratic limb darkening coefficients (71 and 
72); the newly-derived mid-transit times were then used 
to determine an improved orbital period (see Table [1]). 
We also searched for periodicities in the deviations of 
the mid-transit times from our best-fit ephemeris with 
the Lomb-Scargle algorithm, but found no evidence for 
significant transit timing variations. 

The new ephemerides were used to construct a phase- 
folded transit light curve; we then fit a transit model 



to the composite light curve with Levenberg-Marquardt 
least-squares, including a binning of the model to match 
the sampling of the data ()Kippin g 2010). Uncertainties 
calcu lated using the "prayer-bead" method (jGillon et al.1 
2007). The orbital parameters derived from this fit 
are listed in Table [TJ Our linear limb-darkening coef- 
fic ient is consistent with m odel atmosphere predictions 
bv lClaret fc Bloemenl (|2011[ ). but our required quadratic 
limb-darkening coefficient is moderately different than 
their predicted values. 

3.2. Stellar Gravity Darkening 

The residuals of the composite transit light curve of 
HAT-P-7b (see Figure [1]) show an apparent brightening 
of the system just prior to mid-transit. After fitting for 
the transit model, we found x 2 ~ 720. This is greater 
than the expected value for the 642 photometric mea- 
surements with seven fitting parameters {Xideai = 635), 
giving a p- value of 1%. Thus the anomaly in the transit 
light curve is of marginal statistical significance to verify 
with certainty that it has an astrophysical origin. How- 
ever, we will hypothesize here that it is real and test can- 
didate astrophysical phenomena that can cause a similar 
light curve anomaly. 

Transient stellar phenomena such as magnetically- 
driven star spots can be ruled out as the cause of 
the anomaly, since star spots only persist for time 
scales much shorter than this set of observations, which 
spanned two years. To ensure that the anomaly is not 
a few star-spot crossings that became diluted in the av- 
erage, we fit model transit light curves to four subsets 
of the entire observation, and the anomaly was persis- 
tent in each subset. We therefore require a periodic 
stella r dimmi n g asso ciated with the planet's orbital pe- 
riod. I Barnes! (|2009| ) suggests that gravity-darkening of 
rapidly rotating stars can perturb transit light curves; 
however this effect is only significant for stars that rotate 
much more rapidly tha n HAT-P-7, which has v sin i = 3.8 
km/s (Pal et al. 2008) and an ex pected rotation velocity 
v ~ 15 km/s (jWinn et al.1 12009). Furthermore, gravity 
darkening due to rotation affects the light curve on all 
times between the second and third contact, whereas the 
anomaly that we detect is highly-localized in phase angle. 
One possible explanation for such a localized anomaly is 
planet-induced gravity darkening, in which the tidal dis- 
tortion of the star due to the planet causes a decrease 
in the surface gravity near the sub-planet region, which 
would result in a corresponding decrease in stellar sur- 
face brightness (J12). Here, we construct a simple model 
to characterize the plausibility of gravity darkening as an 
explanation for the anomaly in the residuals. 

The dark spot is assumed to follow a Gaussian intensity 
profile with its centroid located in the orbital plane of the 
planet, displaced from the sub-planet point on the stellar 
surface by some phase lag, S , as illustrated in Figure [5J 
We calculate the brightness of the system as the planet 
transits and the dark spot traverses the surface of the 
star. The total flux as a function of orbital phase is 
then compared to a similar model with no dark spot to 
demonstrate the difference between our model and the 
expected transit light curve. The difference predicts the 
amplitude, phase lag and width (a s ) of the anomaly in 
the transit light curve residuals. 

We fit our spot model to the residuals of the tran- 
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TABLE 1 
HAT-P-7 System Parameters 



Parameter 



Definition 



Kepler Photometry 



P (days) Period 
T c (BJD TDB) Epoch 

Inclination 



i (deg) 
Rp/ R* 

7i 
72 

a/R„ 

D (ppm) 
T B (K) 



Ratio of planetary 
radius to stellar radius 

Limb-darkening, linear 

Limb-darkening, quadratic 

Ratio of semimajor axis 
to stellar radius 

Secondary eclipse depth 

Brightness temperature 

Phase at mid-eclipse 



2.204737 ± 0.000017 
2454954.357463 ± 0.000005 
83.111 ±0.030 

0.07759 ± 0.00003 

0.3525 ± 0.0066 
0.168 ±0.010 

4.1502 ± 0.0039 

69.1 ±3.8 
2733 ± 21 
0.500051 ± 0.000006 



Ephcmerides (i.e., P and T c ) derived from Levenburg-Marquardt least-squares fits to each transit, and orbital parameters derived by the 
Levenberg-Marquardt least-squares fit to the composite transit light curve with one-minute binning. The uncertainties of the parameters 
were calculated using the prayer-bead method. The linear and quadratic limb-darkening coefficients (71 and 72) are defined as in the 
IMandel fc Alio] 112003) formalism. 

model. We described the emission spectrum of HAT-P-7 
with that of a blackbody with a n effective temp erature 
of T* = 6350 K (as measured bv lPal et all [2001 . since 
conventional stellar atmospheric models do not simulate 
the differences between distinct atmospheres with tem- 
perature differences on the order of tenths of Kelvin. We 
found a "best-fit" spot model with phase lag of the Gaus- 
sian centroid from the sub-planet point on the stellar sur- 
face 4> s = 0.056 [phase units], with amplitude -106 ppm 
and radial spread <r s = 0.037 [phase units]. We found 
that the most consistent temperature difference between 
the darkened spot and the mean stellar surface, calcu- 




Dark Spot 



111 

To Observer 



lated from this amplitude, is T< 



spot 



T* = -0.183 K, 



Ingress 



Egress 



Fig. 2. — The system geometry is depicted above, with the planet 
represented by the small dark circle, and the star represented by the 
large grey circle, with a gravity-darkened spot on its surface. The 
transit of the planet across the dark spot occurs after ingress and 
before mid-transit, so long as the phase lag between the sub-planet 
point P and the centroid of the spot, rf> B , is small and nonzero. The 
diagram depicts the beginning of the transit across the dark spot, 
viewed at an angle such that the orbital plane is in the plane of the 
page. Note that the plane containing the planet and the observer 
is tilted by 90 — i = 6.9 degrees to the plane of the page. The 
schematic is not drawn to scale. 

sit light curve (Figure [1} to determine the amplitude 
and phase lag of the darkened spot, then we convert 
the amplitude of the anomaly in ppm to Kelvins us- 
ing a blackbody model for the star. The model which 
minimizes the % 2 computed between that model and our 
composite transit light curve is adopted as the "best-fit" 



consistent with the anomaly of "a few 0.1 K" predicted 
by J12. 

Using our derived orbital inclination and a/i?*, we find 
that as the planet transits the spot, the location of the 
sub-observer point (the projection of the transit onto the 
star) is 29.6 degrees above the line of sight to the center 
of the star. If we assume that the path of the center of the 
spot has the same inclination as the planetary orbit (6.9 
degrees from the line-of-sight to the star), the fact that 
the planet transits the spot at all means that the spot 
must have a half-width of > 22.7 degrees in the direction 
orthogonal to the plane of the planet's orbit. We found 
<7 S , the 1-sigma width of the spot deduced from the fit to 
the transit across the spot in the direction of the orbit of 
the planet, to be a s — 0.032, or 11.5 degrees. Assuming 
that the dark spot is real, it is evidently more extended 
in the direction orthogonal to the plane of the planet's 
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orbit. 



3.3. Secondary Eclipse 



As mentioned above, several previous studies have 
measured the mean depth for the secondary eclipse of 
HAT-P-7b, but t hese results are all sign ificantly different 
from each other. iBorucki et al.l (2009) fo und a n eclipse 
depth of 130±11 ppm. while lWelsh et all (f20Toh and J12 
found secondary eclipse depths of 85.8 ppm and 61 ± 3 
ppm, respectively. Each study used a different number 
of observations, as well as different methods to compen- 
sate for the planetary phase variations and the stellar 
ellipsoidal variations. With the goal of removing model 
uncertainties, we masked the eclipse and fit model phase 
curves to the light curve to account for the significant 
phase variation effects, Doppler beami ng, planetary re- 
flecti on and ellipsoidal variation curves (jFaigler fc Mazehl 
l2011f l. from phase 0.3 to 0.7. We then divided by this 
baseline before fitting to the eclipse. We hereafter chose 
to present results from the PDCSAP fluxes, which are 
consistent with the results f r om th e SAP fluxes. 

We fit a iMandel fc Ago] (|2002f ) model eclipse light 
curve to each eclipse individually using Levenberg- 
Marquardt least-squares. The orbital parameters found 
by the previous fit to the composite transit light curve 
(see Figure [T]) were used to fix the duration of the eclipse 
and the mid-eclipse times; uncertainties for each depth 
measurement were estimated by the prayer-bead method. 
We then derived the mean eclipse depth by computing a 
least-squares fit to the eclipse depth over all phases in- 
cluded in our analysis of Quarters 1-9, yielding a depth 
of 69.1 ± 3.8 ppm (see Figure [3]). The mid-eclipse time 
derived from the fit to the composite secondary eclipse 
light curve occurs virtually exactly at the expected value 
(phase = 0.5), after taking light travel time across the 
orbit into account. This va lue is significantly low er than 
the eclipse depth found by IBorucki et al.l (|2009| ). due to 
both the longer observational baseline and the inclusion 
of ellipsoidal va riations in our mode l; fitting only the 
data analyzed bv lBorucki et al.l ()2009h . we find an eclipse 
depth of 105 ± 27 ppm. 

The individual eclipse depths exhibit scatter of 32 ppm 
over the full data set, assuming a constant expected 
value. This scatter is higher than can be accounted 
for using the Kepler -provided uncertainties (designated 
"SAP .ERR" and "PDCSAP_ERR"). In order to ensure 
that the scatter we measured does not include intrinsic 
variability in the eclipse depth, we repeated the analysis 
to search for an eclipse at phase=0.3, where no eclipse 
is expected. Our analysis is consistent with no eclipse 
at phase=0.3 with 32 ppm rms scatter, which confirms 
that the scatter between measurements is indeed non- 
astrophysical. This complete analysis was repeated using 
both the SAP and PDCSAP fluxes and similar measure- 
ments were obtained for the eclipse depth and scatter. 

Scaling the geometric albedo result from C10 with our 
smaller measurement of the eclipse depth, we find a geo- 
metric albedo of < 0.03. This low optical albedo of HAT- 
P-7b is consistent with observations of other highly irra- 
diat ed hot Jupiters, suc h a s TrES-2b with albedo ~ 0.01 
(see lBarclav et al.l2012landlKipping fc Spiegelll2011l ) and 
HP 209458 (|Rowe et all 12008ft. Theoretical predictions 
also call for low albedos for hot Jupiters, such as the 
atmospheric models by C10. 



To compute the brightness temperature of the planet, 
a Kurucz model star with temperature T„ = 6350 K is 
adopted for HAT-P-7. If we assume that all of the flux 
from the atmosphere of the planet is thermal emission, 
the ratio of the blackbody flux of the planet to the Ku- 
rucz flux of the star integrated over the Kepler bandpass 
gives the secondary eclipse depth. These considerations 
yield a brightness temperature in the Kepler bandpass of 
Tb = 2733 ± 21 K, where these uncertainties propagate 
from the rms scatter of the eclipse depth (see Figure [3]) . 

4. DISCUSSION 

To further test the hypothesis that the anomaly is 
caused by a gravity-darkened spot on the stellar surface, 
we compare the physical characteristics of the dark spot 
to our knowledge of stellar atmospheres. The dark spot 
is offset in phase from the sub-planet point on the stel- 
lar surface by phase=-0.059. In other words, about 3.1 
hours elapse at the sub-observer point on the stellar sur- 
face between the passing of the planet overhead and the 
passing of the centroid of the dark spot. 

The sound speed in the photosphere sets the limit on 
how quickly the stellar atmosphere can react to the grav- 
itational perturbation of the passing planet, and thus the 
minimum elapsed time expected between the passing of 
the planet and the passing of the dark spot through the 
sub-observer point is of order 



At 



H 



(1) 



where H is the scale height of the atmosphere, and c s is 
the sound speed. Using a perfect gas equation of state, 
we find 



At 



' 3k B T 
5fim H g 2 



(2) 



The temperature (T) and surface gravity (g) at the pho- 
tosphere, where the dark spot manifests, can be inferred 
from simulations of stellar atmospheres to calculate At. 
Robert Kurucz has calculated the pressure, temperature, 
density and sound speed at many depths in simulated 
stellar atmospherefl We chose the nearest analogue in 
the grids of simulated stars to HAT-P-7 with surface tcm- 
perature TeJt = 6250 K and surface gravity logg = 4 
(|Pal et a"H 12008ft . We compute the temperature of the 
photosphere, which we define as the layer of the atmo- 
sphere where the optical depth r ~ 2/3. Our observed 
response time of 1.06 x 10 4 s, measured by the phase lag 
of the dark spot, is essentially equivalent to the Eq. [2] 
prediction (At ~ 1.02 x 10 4 s), supporting our hypoth- 
esis that the perturbation seen in Figure [1] is due to the 
time-dependent gravity-darkened response of the stellar 
atmosphere to the passage of the planet, as predicted by 
J12. 

5. SUMMARY 

We have presented new analysis of eight quarters of 
Kepler photometry of HAT-P-7b. We refine the orbital 
parameters of the system and a smaller secondary eclipse 
depth than the earliest published results. These findings 

3 Models by Kurucz available on his website: 
http : / /kurucz . harvard . edu/gr ids . html 
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Secondary Eclipse Depths (Kepler Ql-9} 
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Composite Secondary Eclipse (Kepler Ql-9) 
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Fig. 3. — Upper: Eclipse depth measurements for each phase 
analyzed, with best-fit depth = 69.1 ± 3.8 ppm (horizontal line). 
This depth is consistent with predominantly thermal emission, with 
nearly no reflected light component. Lower: Phase-folded eclipse 
light curve for Kepler Quarters 1-9. The photometric measure- 
ments are averaged in ~ 1.5 minute bins (cir cles). Error-bars are 
omitted for clarity. The curve represents the Mandcl & Agol (2002) 
model light curve for a secondary eclipse of a uniformly luminous 
occulted planet. 

are consistent with other analyses of subsets of the ob- 
servations that we considered. The secondary eclipse of 
HAT-P-7b indicates a geometric albedo of < 0.04, quali- 
tatively consistent with the properties of the often com- 
pared, highly-irradiated hot Jupiter TrES-2b. We find 
no evidence for variations in the eclipse depth. 



The composite transit light curve of HAT-P-7b has 
a ~ 10 ppm anomaly that is consistent with both the 
amplitude and phase lag of a planet-induced gravity- 
darkened spot on the stellar surface, though additional 
transit observations are required to confirm it. The sub- 
planet region of the star experiences a decrease in surface 
gravity, which results in a corresponding decrease in tem- 
perature and brightness that is observable in the light 
curve. The ongoing Kepler mission will continue to ob- 
serve HAT- P- 7b, and future analyses of the phase-folded 
transit light curve will be able to conclusively verify the 
astrophysical nature of the anomaly. 
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